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ABSTRACT 

Observations of stellar prominences on young rapidly rotating stars provide unique probes 
of their magnetic fields out to many stellar radii. We compare two independently obtained 
datasets of the K3 dwarf Speedy Mic (BO Mic, HD 197890) using the Anglo- Australian Tele- 
scope (AAT) and the ESO Very Large Telescope (VLT). Taken more than a fortnight apart 
\Q ' they provide the first insight into the evolution of the prominence system on such a young 

rapidly rotating star. The largest prominences observed transiting the stellar disc are found at 
very similar rotational phases between the epochs. This suggests that the magnetic structures 
supporting the prominences retain their identity on a two to three week timescale. By taking 
advantage of the high signal-to-noise and large wavelength range of the VLT observations 
we identify prominences as transient absorption features in all lines of the Hydrogen Balmer 
series down to Hio. We use the ratios of the prominence EWs in these lines to determine 
their column densities in the first excited state of hydrogen. We determine the optical depths, 
Q-T finding prominences to be rather optically thick (r « 20) in the H Q line. The total hydrogen 

column density and thus the prominence masses are determined via observations of the Call 
■ H&K lines. We find typical masses for four of the largest prominences to be in the range 

^ \ 0.5 — 2.3 xlO 14 kg, slightly larger than giant solar prominence masses. Rotationally modu- 

Ci ■ lated emission is seen outside of the H Q line. These loops of emission are shown to be caused 

by prominences seen off the stellar disc. We find that all of the large emission loops can be 
associated with prominences we see transiting the stellar disc. This combined with the fact 
. that many prominences appear to eclipse the off-disc emission of others, strongly suggests 

$_i ' that the prominence system is highly flattened and likely confined to low stellar latitudes. 
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1 INTRODUCTION 

Solar prominences are either observed against the surface of the 
Sun as dark absorption structures known as 'filaments' or off the 
solar limb as bright emission loops. Traditionally to observe such 
features a narrow-band filter is used centred on 6563A. This corre- 
sponds to the first and strongest transition of the hydrogen Balmer 
series, hydrogen alpha (H Q ). 

Stellar p rominences were first reported by 
ICollier Cameron & Robinsor] ll989al) as rapid transient ab- 
sorption features passing through the Doppler broadened H a line 
of the rapidly rotating K0 dwarf AB Doradus (AB Dor). From 
subsequent observations these clouds of material were shown 
to be in presumably magnetically enforced co-rotation with the 
stellar surface and so dubbed "prominences". Since their discovery 
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prominences have now been reporte d on around a dozen rapidly 
rotating late-type stars (see ICameron et alJ 120031 for a recent 
review). These observations have been primarily restricted to the 
strong Balmer lines, often only (H a ) and occasionally hydrogen 
beta (H^), 

The AB Dor prominence system remains the 
most studied to d ate. Pr ominence s have consistently 
been reported {C ollier Cameron & Robinso nl Il989at 
[Collier Cameron & Robinsonl |l989b|: |Donati & Collier Ca meron 
Il997t ICollier Cameron et alj|l999t iDonati et alJll999l) to lie aTa 
range of distances from the stellar rotation axis from 2 to 8 R*, 
with a concentration at the co-rotation radius (2.7 R* on AB Dor). 

iDunstone et alJ {2006) (hereafter Paper I) presented observa- 
tions of the prominence system of another rapidly rotating star, 
Speedy Mic (BO Mic, HD 197890). Speedy Mic is a K3V dwarf 
with a 0.38 d rotational period, an equatorial rotation velocity of 
vsini ~ 132 km s _1 and is seen at a high inclinat i on. Th e stel- 
lar inclination was found to be 70 ± 5° by iBarnesI 120051) . how- 
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ever at such high inclinations Doppler imaging provides poor con- 
straints. IWolter et alJ ^2005) also found an inclination of 70° , yet 
they concede that varying the inclination by up to ±20° has lit- 
tle influence on the resulting Doppler maps. Recent studies of X- 
ray variability have shown Speedy Mic to be one o f the m ost ac- 
tive stars in the solar neighbourhood. ISingh etail Q999) found 
log(Lx I Lboi) — —3 placing it in the saturated regime. In Pa- 
per I we found a densely packed co-rotating prominence system, 
with heights of prominences grouped around 3 R» from the stel- 
lar rotation axis. This meant that they were not concentrated at the 
co-rotation radius (1.95 R* on Speedy Mic) but in fact at twice 
the height above the stellar surface. Clearly such observations of 
co-rotating prominences so far above the co-rotation radius require 
them to be supported by closed magnetic loops. Then prominences 
should be direct tracers of the stellar magnetic field on the largest 
scales, out to many stellar radii in the corona. 

X-ray and Zeeman-Doppler Imaging (ZDI) observations are at 
odds with the hypothesis that large magnetic loops extend high into 
the co rona of such rapidly rotating stars. Recently Hussain et al. 
(2005) showed that the closed field containing compact hot X- 
ray emitting material on AB Dor was close to the stellar surface 
and unlikely to extend beyond a radius of 1 .75 R» . From ZDI ob- 
servations of the surface magnetic field of AB Dor the structure 
of t he X-ray emission was d etermined assuming a simple heating 
law Ijardine. Collier Cameron & Donatil2002tljardine et all2002l) . 
They found that the X-ray emission must be confined fairly close to 
the stellar surface. As Speedy Mic is a faster rotator than AB Dor 
one would expect the coronal extent to be smaller. These obser- 
vations are in apparent conflict with the requirement for magnetic 
support of prominences at heights of several stellar radii. 

A model for prominence support which does not require 
closed magnetic field lines to extend to seve ral stella r radii abov e 
the surface was put forward bv ljardine & van Ballegooiie nl <2005h . 
Instead they propose that prominences are supported by the re- 
connection of oppositely directed wind-bearing field lines over 
coronal helmet streamers. This would result in long thin, highly 
curved, loops extending out to large radii. These could then fill up 
with gas that was originally part of the win d, forming the observed 
prominences. However while Ijardine & va n Ballegooiien 1 2005) 
show that such magnetic loops can achieve equilibrium they do not 
attempt to model the prominence structures themselves. In order to 
test such models we need more information on the physical proper- 
ties, evolutionary scenarios and lifetimes of prominence systems. 

Solar prominences are also be observed in the H a line as emis- 
sion features seen off the solar limb. Photons at this wavelength 
originating from the surface of the Sun can be absorbed by a hy- 
drogen atom in the prominence, thus exciting its electron from 
the first (n=2) to the second (n=3) excited state. Providing such 
an atom does not undergo a collision, the electron will de-excite 
either to the ground state (n=l) via the emission of a Lyman beta 
(Ly^) photon or, almost equally likely, back to the first excited state 
via the release of another H Q photon. Lya photons are re-absorbed 
almost immediately, creating new H a photons. As this emission is 
isotropic most photons originating from the part of the solar surface 
occulted by a prominence are 'scattered' out of our line-of-sight, 
creating an absorption feature. Similarly, photons can be scattered 
into our line-of-sight by a prominence seen off the solar limb, cre- 
ating an emission feature. Observing a stellar prominence as both 
an absorption and an emission feature could potentially provide im- 
portant insights regarding its physical and geometrical properties. 

The initial motivation for this study arose from a fortu- 
itous coincidence. Two weeks after the AAT Speedy Mic obser- 



vations used in Paper I, a separate team performed very simi- 
lar high-resolution spectroscopic observations of Speedy Mic with 
the VLT. Two papers using this da t aset h ave been pu blished by 
the or i ginal observers. IWolter et alJ |2005I) and IWolter &■ Schmifl 
l2005h1woTter et alJ (2005) provides a full account of the obser- 
vations and then goes on to use the red-arm for Doppler imaging. 
IWolter & SchmittN2005h focus on the potential of the Call H & K 
lines for future Doppler imaging of the stellar chromosphere. 

Here we compare the two sets of observations in [J2| then in Sj5| 
we use the VLT observations to obtain prominence physical prop- 
erties. In S|4| we report the observation of loops of emission and 
show that these are caused by prominences seen off the stellar disc. 
We discuss our results in the context of the coronal structure of 
Speedy Mic in iJ5] Finally we present our conclusions in f|6| 



2 OBSERVATIONS 

The Speedy Mic observations were obtained from the ESO Science 
Archive Facility. They consist of two nights (2002 August 02 and 
07) of data using the 8-m Very Large Telescope (VLT) UT 2. The 
spectra cover a large wavelength range by using the UVES spec- 
trograph in dichroic mode, 3260A to 4450A (the "blue-arm") and 
4760A to 6840A (the "red-arm"). A l"slit width gave a spectral 
resolution of A/AA « 40 000. Speedy Mic was observed contin- 
uously on both nights yielding complete and uninterrupted phase 
coverage. 

We re-extracted the data using ECHOMOP, the echelle re- 
duction package developed by Mills 1 1994) and its implementation 
of the optimal extraction algorithm developed by iHornd ll98d) . 
During the data reducti on we encountered the same problem re- 
ported bv lWoltereta?] <2005l) of variations in the MIT-CCD red- 
arm continuum, presumably due to the CCD manufacturing pro- 
cess. In order to correct for the introduced small scale structure we 
created a "balance frame" that included these variations. 

In total there are 125 spectra on August 02, each with an ex- 
posure time of 200 s. The readout time of 70 s makes the cadence 
of exposures 270 s, therefore 122 exposures were required to cover 
the 0.38 day period of Speedy Mic. The situation on August 07 is 
slightly more complicated with the exposure times being adjusted 
between 120 s and 200 s to compensate for variable sky conditions. 
Four spectra had to be abandoned during reduction due to over- 
exposure, yielding a total of 141 usable spectra. 

After each spectrum was extracted we normalised the spec- 
tral orders by fitting the continuum using splines. Due to the dif- 
ficulty in locating the continuum level in the blue orders we used 
a two-dimensional fit that constrained the local shape of the con- 
tinuum. We discuss the limitations of this technique in J3| The re- 
sulting signal-to-noise (S/N) of individual spectra were very high. 
For example near the centre of the order containing the H a line the 
unbinned S/N ranges between 350 and 600 per 0.05A wavelength 
element throughout the night. 



2.1 A fortnight of evolution 

Thanks to the fortuitous coincidence of two different Doppler 
imaging groups (using different telescopes) both observing 
Speedy Mic within two weeks of each other we now have an excit- 
ing and unique opportunity. The two to three week time gap is not 
normally an easy one to achieve with high-resolution spectroscopy. 
Practical considerations like instrument availability on any given 
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Figure 1. Raw H Q timeseries spectra for 2002 July 18, August 02 & 07, with phase plotted against velocity. Dashed black lines show the v sini limits and the 
central dashed white line shows the radial velocity of Speedy Mic. The grey-scale runs from black at 0.85 times the local continuum level, to white at 1.05 
times continuum. The large prominences discussed in the text are labelled. 



telescope normally limit the time span of observations at any given 
epoch to mere days. 

In FigQwe have plotted the H a trailed spectrum (timeseries) 
of July 18 from Paper I and compare it to August 02 and 07. Note 
that for comparison with this work the July 1 8 timeseries has been 
re-phased taking the epoch of phase zero to be the first observation 
on August 02. There are fifteen days between July 18 and August 
02, that equates to approximately 40 stellar rotations and there are 
twenty days or 53 rotations between July 18 and August 07. Despite 
the differences in the quality of the data between the two epochs it 
is clear to see that qualitatively the timeseries appear similar. 

A visual examination of FigQshows that August 02 does in- 
deed share features of both July 18 and August 07. We can see that 
one half (the 'dark' side) of the Speedy Mic rotation still has in- 
creased absorption due to large prominences. In the other half (the 
'light' side) of the stellar rotation cycle the surface is more easily 
seen and we observe weaker prominences embedded in and around 
bright surface 'plage' features. Identifying and comparing individ- 
ual prominences is difficult as there has been considerable evolu- 
tion between the July and August epochs. Very few of the weaker 
prominences observed on the light side can be re-identified. The 
large dark region at <f> « 0.35 appears stronger and more compact 
in phase in the August observations. The strong group of promi- 
nences at <f> ~ 0.15 seen on both August epochs has no counterpart 
in the July observations. 

The only prominence that can be compared with any confi- 
dence is prominence A at <j> « 0.55. This prominence was reported 
in Paper I (again labelled 'A') as being present on all three nights 
(July 18, 19 and 23). On August 02 we find a prominence at a very 
similar phase. Yet on August 07 there is no evidence of the promi- 



nence - only a bright surface region. So can this prominence seen 
on August 02 really be the same structure that was observed fifteen 
days earlier? The main properties of a prominence that we can mea- 
sure at both epochs are its location in the rotational phase, the pro- 
jected distance from the Speedy Mic rotation axis (height) of the 
prominence and the maximum fractional depth below the contin- 
uum. The first two of these are found from the Gaussian matched- 
filter tracking technique described in detail in Paper I. The differ- 
ence in the two phases of observation between July 18 and Au- 
gust 02 is A<j) — 0.007, which is the same as the phase separation 
between consecutive exposures. The height of prominence A was 
quoted as 2.79 ± 0.11 for July 18 and is found to be 2.99 ± 0.07 
on August 02. The maximum fractional depth relative to the con- 
tinuum when the prominence is transiting the disc center is « 0.87 
on July 18 and « 0.88 on August 02. 

The similarity of the properties of this prominence at both 
epochs lead us to conclude that this is almost certainly the same 
structure. As we shall discuss in [J5| however, we are not suggest- 
ing that this prominence necessarily contains the same material at 
both epochs. 



2.2 Spectral lines showing prominence signatures 

After each spectral order was extracted we stacked the exposures to 
create a timeseries for each order. We then collapsed this timeseries 
to create a mean spectrum and subtracted it from each individual 
spectrum. This left us with only time varying features, i.e. star- 
spots, surface plage regions and prominences. 

We performed a visual inspection of all orders for evidence 
of prominence absorption features moving rapidly through the line 
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profile, just as we observe in the H Q line. This dataset is ideal for 
an unbiased search, as it is of exceptional S/N and the large wave- 
length range gives us plenty of lines to examine. As was expected, 
the lines of the Hydrogen Balmer series all show evidence of promi- 
nence material. Starting with the strong K a line the same charac- 
teristic prominence pattern is repeated, with different strengths, in 
each line. It is possible to visually identify this pattern all the way 
down to H14 at A = 372lA. Prominences in the later hydrogen 
lines are so optically thin that only the very strongest prominence 
is identifiable. Furthermore the star-spot residuals from the many 
photospheric metal lines present at these blue wavelengths have 
strengths of the same order of magnitude as the prominence signa- 
tures. Therefore for the rest of this paper we restrict our analysis to 
lines down to Hio at A = 3835A, all of which are shown in Fig|2| 
The singly ionised calcium (Call) H & K lines at A = 3968A 
and A = 3933A have strong emission cores and show strong 
prominence absorption signatures. These lines are discussed at 
length in 33.21 Finally, we find very weak prominence signatures 
in the sodium doublet at A = 5890A. Unfortunately these are 
poorly placed at the very start of an order; this combined with their 
strengths make them unsuitable for further analysis. These were the 
only lines that we can say with confidence show any prominence 
material. 



3 PROMINENCE MASSES 

3.1 Column density of hydrogen in the n=2 level 

Due to the large wavelength range and the high S/N of these obser- 
vations it is possible to compare the strengths of individual promi- 
nences in different hydrogen Balmer lines. In order for such a com- 
parison to be possible each line has been re-binned to the same 
velocity scale, 3.1 km s _1 per pixel. 

The trailed spectra of all eight Balmer lines being considered 
are shown in Fig. [2] for August 02. Some of the lines have small 
problems with their reduction. Hp clearly has four bad spectra, 
which arise due to a bad column in the earlier part of this spec- 
tral order. When ECHOMOP traces this order it rejects and inter- 
polates over these bad pixels, however this procedure was unsuc- 
cessful for these four spectra and a good trace was not obtained. 
Consequently they are not used in further analysis. Hio shows a 
vertical line which is due to another chip defect. Finally it should 
be noted that the H E line is superimposed with the Call H line, due 
to the strong prominences seen in Call we do not use H e further in 
this analysis. 

3.1.1 Calculating the curve of growth 

The ratio of the equivalent widths (EWs), or strengths, of individual 
prominences in different hydrogen lines can be used to calculate the 
column density of hydrogen atoms in the n — 2 level. 

In order to do this we need to first calculate the absorption 
cross-section of a hydrogen atom at the centre of each Balmer line. 
This is given by: 



where g u and gi are the statistical weights of the upper and lower 
levels of the transition, c is the speed of light and v is frequency. 
A u i is the Einstein A-coefficient and gives the probability of spon- 
taneous decay from the upper to the lower level. $ is the line profile 



for which we use a Voigt profile approximation. The Lorentz pa- 
rameter was taken to be natural (radiative) dampening. This is jus- 
tified because of the scattering mechanism that causes us to observe 
the prominences in absorption whilst transiting the stellar disc. This 
tells us that the time between collisions within the prominence is 
longer than the spontaneous de-excitation time and so collisional 
broadening is unimportant. We assume the dominant broadening 
comes from two sources. Firstly the thermal Doppler velocity of 
hydrogen at a prominence temperature of 10,000 K which is 12.9 
km s _1 and secondly a random turbulent velocity. 

Solar prominences exhibit random turbulent motions at least 
as great as a few km s _1 (Engvold 1998). The amount of turbu- 
lence we include in the line profile will have a greater effect on 
the Call H&K observations (discussed in 33,21 than on those of 
the hydrogen Balmer lines. This is because the calcium atom is 
heavier than the hydrogen atom and so has a much smaller thermal 
Doppler velocity of 2.0 km s _1 . We can place an upper limit on the 
possible turbulent velocity in Speedy Mic prominences by measur- 
ing the widths of the smallest prominences in the Call H&K lines. 
A width of approximately 18 km s _1 is found for the smallest ab- 
sorption features as they cross the centre of the stellar disc. If we 
assume that these small prominences are unresolved then the mea- 
sured width should be due to turbulence, but we also need to con- 
sider the rotational smearing that occurs during a 200 s exposure. 
The average prominence height found in Paper I of 3 R* corre- 
sponds to an orbital velocity of 400 km s _1 . The Speedy Mic ra- 
dius was calculated in Paper I from a simple dynamical argument 
to be 1.06 ± 0.04 R.0, considerably more bloated than its K3V 
main-sequence radius. Therefore during a single exposure such a 
prominence travels approximately 11% of the Speedy Mic radius 
which in velocity space is 14 km s" 1 , the majority of the observed 
width. If we assume that these velocities add in quadrature then we 
arrive at an upper limit of 11 km s _1 for the turbulent velocity. The 
validity of these assumptions are discussed further in 35.21 Now we 
assume a turbulent velocity of 5 km s _1 for the following results. 

The optical depth is found by multiplying the absorption 
cross-section per atom by the column density (N2 - the total num- 
ber density of hydrogen atoms in the n — 2 level per square meter): 

T„ = cr„ N2 

Therefore the EW is found for each line by integrating over 
the line profile: 



EW = 




By considering a range of different values for the column den- 
sity of atoms in the n = 2 level we obtain a curve of growth for 
each line. This describes how the EW of the line changes as we 
increase the density of absorbers. At low densities (optically thin) 
the curve of growth is a linear increase in the EW proportional to 
the probability of the transition. As we increase the column den- 
sity eventually each line starts to saturate in its core leading to a 
slowing in the growth of the EW. When the Doppler core is totally 
saturated (optically thick) then further growth of the EW occurs at 
a rate determined by the opacity in the Lorentzian line wings. 

Once we have the EW of each line as a function of column 
density it is then possible to express each line's EW as a fraction 
of EW(H a ). These ratios are shown as the solid lines in the left 
panel of Fig|3| This shows that H Q starts to saturate first. Its growth 
of EW slows and the relative strength of the later lines increases. 
One by one the other lines start to saturate also. Therefore by mea- 
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Figure 2. The first eight lines of the hydrogen Balmer series for August 02 are shown. An average profile has been subtracted from each timeseries (see text). 
The contrast is such that black is 0.85 and white is 1.05 times the local continuum level. Dashed vertical lines show the v sin i limits. Note H f is blended with 
the CallH line and Hp has four bad spectra (see text). 
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suring the ratio of the observed EWs from the timeseries it should 
be possible to find a unique value for the column density of each 
prominence. 

3.1.2 Measuring the EW ratios 

As mentioned earlier it is unfortunately impossible to use H e in 
this analysis due to it being blended with the Call H line. Even 
those parts of the line not directly superimposed on the calcium 
line could be subject to off-disc emission from it (as discussed in 
S|4}. We also had to abandon measuring the Hg line due to anoma- 
lous results caused by a very strong neighbouring metal line. The 
star-spot residuals in this line appeared strong enough to affect our 
prominence EWs. Thus the analysis is now restricted to: H a , H^, 
H 7 , Hi, Hg and Hio. 

To measure the EW of the prominence absorption signatures 
it is first necessary to remove the 'time-invariant' background stel- 
lar spectrum. For Speedy Mic this is a rather more difficult task 
that one might first anticipate, the problem being that much of the 
timeseries is in strong absorption due to the prominence material. 
Other parts of the timeseries, however, show strong surface emis- 
sion features (plage regions). The nature of the background sur- 
face spectrum for the regions covered by large prominences is un- 
known from H Q alone. It could well be that these phases are in even 
stronger emission. By examining the later Balmer lines in Fig. [2] 
where even the strongest prominences become optically thin, we 
can see that this is probably not the case. It is also important to re- 
member that the bright surface plage regions contribute to the flux 
that the prominences receive and so should not simply be ignored. 
We finally chose to represent the background spectrum by taking 
the mean of the 'light' half of the timeseries which is devoid of 
strong prominences. Specifically the phases 0.65 to 1.0 were used 
to create the mean profile that was then subtracted from the entire 
timeseries. Initially we tried a much smaller range of phases which 
apparently had no prominences. However this produced an asym- 
metric mean profile and upon further examination a bright surface 
feature was present. Using this mean profile we obtained column 
densities that differed by a factor of 1.25 to 2.5 from those reported 
below. This is similar to the level of uncertainty we find on the 
prominence column densities. 

We use the high S/N H Q timeseries to identify which pixels 
in the timeseries are to be used for the analysis. In practice this 
was done by selecting a box around a single prominence feature 
and selecting all the negative pixels within that box. The EWs are 
then simply the ratio of the sum of the corresponding pixels in each 
line to that of H a . Associated uncertainties on these ratios are also 
derived. 

The points where the observed EW ratios intersect the theoret- 
ical ratio of the curves of growth are shown for each line in Fig. [5] 
The uncertainties in measured EWs are then translated into uncer- 
tainties in column density using the model curves of growth ratios. 
Thus we can take an inverse variance weighted average of these 
points to obtain a single estimate for the column density in the n=2 
level of hydrogen. The final uncertainty on the column density of 
each prominence is estimated from the standard deviation of the 
different line ratios. 

Obviously the high S/N of the stronger Balmer lines like Hp 
allow us to measure the ratio of their EWs to H a to a much greater 
precision than some of the weaker, later, Balmer lines. However 
because it turns out that the cloud is quite optically thick in Hp 
it means that there is a larger uncertainty on the column density 
once the measurement uncertainty has been converted to an un- 



certainty in the column density. Therefore H 7 and Hg have most 
weight when calculating the optimal column density. 

This technique only works for the very strong prominence fea- 
tures. As already discussed there are many weaker prominences 
that cross the star's brighter hemisphere (phases 0.6 - 1.0), most of 
which can only be detected in H a and Hp. We therefore limit our 
analysis to the three strongest prominence signatures on August 02 
(labelled A,B and C in Fig0 and the strongest prominence on Au- 
gust 07 (labelled D in Fig[T). 

The results of this analysis are shown in Table Q We find 
that the prominences have a relatively small spread of column 
densities in the n = 2 level of hydrogen, with the mean being 
log N2 = 18.61 ± 0.05 m -2 . Prominence A has the highest col- 
umn density at log N3 = 18.68 ± 0.19 m -2 . As the uncertainty 
on each individual prominence is greater than the scatter in ob- 
served prominence column densities, we quote a final representa- 
tive column density of log N2 = 18.61 ± 0.2 m~ 2 . Multiplying 
the column densities by the central absorption cross-section in each 
Balmer line we find the optical depths of prominences in each line. 
In Table Q we quote optical depths in H a for each prominence. 
These range between 19 < r < 26 confirming that the promi- 
nences are indeed optically thick at the H a line centre. The optical 
depth drops to unity at around H^ or Hg depending upon the indi- 
vidual prominence. 

For the results shown in Table Q we assumed a turbulent ve- 
locity of 5 km s . As this is considerably smaller than the ther- 
mal Doppler velocity of 12.9 km s _1 turbulence had little effect on 
the resulting column densities. If we increase the turbulent velocity 
to the upper limit of 11 km s _1 the column density increases by 
log A r 2 = 0.1 m~ 2 (or a factor of 1.25). This is smaller than the 
uncertainty quoted on each prominence column density in TableQ 



3.2 Column density of Call 

In order to find the masses of the prominences it is necessary to 
look at resonance lines, that is transitions down to the ground state. 
While we were able to use many Balmer lines to estimate the col- 
umn density in the n=2 level, the number of atoms in the ground 
state (n=l) is so far unknown. The ratio between the populations of 
the n=l to n=2 depends on the excitation of the cloud. Two strong 
lines of singly ionised calcium, Call H&K, are transitions down to 
the ground state. 

We plot the Call H&K timeseries spectra in Fig. [5] lust as 
for the each of the Balmer lines in Fig. [2] we have subtracted a 
mean stellar spectrum obtained from the same range of phases. The 
overall general appearance of the prominence pattern looks quite 
similar to that in the Balmer series. One half of the timeseries 
shows dark prominence absorption signatures crossing the stellar 
disc. The other half shows bright plage active regions on the stellar 
surface. Also visible in the timeseries is a stellar flare at <j) ~ 0.3. 

This Call H &K dataset has already been published by 
Wolter & Schmitt| l2005h in the context of using Doppler imaging 
to map the Speedy Mic chromosphere. They refer to unexplained 
'changes on rapid timescales'. As can be seen from Fig[5]these are 
almost certainly due to the prominence material transiting the stel- 
lar disc. 



3.2.1 Calculating the curve of growth 

We again use the ratio of the measured prominence EWs in order 
to ascertain the column density in Call. The curves of growth are 
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Figure 3. Obtaining column densities in n=2 of hydrogen (left panel) and Call (right panel). Solid lines show theoretical curve of growth EW ratios, as com- 
puted in i|3.1.1l for an example turbulent velocity of 5 km s — 1 . Circles show intersection with observed ratios for prominence A and triangles for prominence 
B. The dashed vertical line is the optimal average column density for prominence A and the dot-dashed line for prominence B. Note H e is absent due to 
blending with Call H line and that Hs is anomalous (see text). 



Table 1. The measured values and derived quantities are display along with their uncertainties. These results assume a turbulent velocity of 5 km s — 1 (see 
text). Note that there is a systematic uncertainty on the Call column density of log A r CaiY=±0.25 (see text). We also report the average and standard deviation 
of the four prominences for various physical properties. 



Prominence 


Fractional 


Optical depths 


Column densities [m 


- 2 ] 


Areas 


Masses 




Depth 


T(H a ) 


r(CaII K) 


log 7V 2 


log N CaII 


log iVi 


log A [m 2 ] 


log M [kg] 


A 


0.12±0.015 


25.7 


21.3 


18.68±0.19 


17.88±0.25 


23.53±0.25 


17.31±0.06 


14.06±0.26 


B 


0.25±0.005 


19.5 


15.8 


18.56±0.21 


17.75±0.25 


23.40±0.25 


17.63±0.03 


14.25±0.25 


C 


0.19±0.01 


22.4 


13.5 


18.62±0.19 


17.68±0.25 


23.33±0.25 


17.51±0.04 


14.06±0.25 


D 


0.24±0.01 


21.9 


11.7 


18.61±0.18 


17.62±0.25 


23.27±0.25 


17.61±0.04 


14.10±0.25 


Averages 




22.4 ±2.6 


15.6 ±4.2 


18.61 ±0.05 


17.73 ±0.11 


23.38 ±0.11 


17.52 ±0.15 


14.12 ±0.09 



calculated in the same manner as for hydrogen. As in ^3 .1.11 the 
turbulent velocity is fixed at 5 km s _1 for the following work. 

The oscillator strength of the Call K line is exactly twice that 
of the Call H line. Therefore when examining optically thin mate- 
rial the EW of Call H will be half that of Call K. As we increase 
the column density and the Call K line starts to saturate, the rate 
of EW growth slows, so the ratio of the strengths of the two lines 
decreases. This is shown in the right-hand panel of Fig. [3] 



3.2.2 Measuring the EW ratios 

In order to measure the EWs both Call lines were binned onto the 
same wavelength scale as the Balmer lines. The same points in the 
timeseries identified from the H a timeseries were used to find the 
EWs. The only difference, however, is that we only used the left 
hand side, the side with negative velocity with respect to the centre 
of both lines. This is because of the H E contamination of the Call 
H line. This extends to just under half of the Call H line profile. To 
investigate what effect this may have on our results we performed 
the same analysis on the hydrogen Balmer lines and compared it 
to our results in i|3. II which used the whole of the timeseries. We 
found that this produced a scatter in the derived column densities 
of up to log N2 = 0.2 m -2 (or a factor of 1.6). However the order 
of column densities remained unchanged, e.g. prominence A was 
still the most optically thick prominence. 




£ I , , , , i i i i i i i i i i i i i i i 

"0 5 10 15 20 

Turbulent velocity (km/s) 

Figure 4. The Call column density is plotted as a function of the turbulent 
velocity for prominence A. The axis on the right of the plot shows the cor- 
responding optical depth of the Call K line. The vertical dashed line at 1 1 
km s — 1 gives an observational upper limit on the turbulence. The dot-dash 
horizontal line shows the optical depth of the H Q line for reference. 



From Fig|3|there appears to be an ambiguity on the Call col- 
umn density due to there being two potential solutions. The results 
using the first (lower column density solution) are shown in Table 
\l\ A similar spread is observed in the column densities of Call as in 
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the Balmer lines with a mean of (log Noall = 17.73±0.11 m~ 2 ). 
Prominence A is also found to have the highest column density at 
log N2 = 17.88 m -2 . We also quote in TableQthe optical depth 
of each prominence in Call K, finding 12 < r < 21, less than 
the optical depths in H a (19 < r < 26). This is actually similar 
to the optical depths in H^. Visual comparison of the line in 
Fig|2|with the Call K line in Fig|5|shows the prominence pattern to 
be similar in strength. This is why we can safely exclude the sec- 
ond solution shown on Fig|3| The second solution gives a column 
density approximately a thousand times higher than that of the first 
solution. Therefore the prominence optical depths in Call K would 
also be a thousand times higher. Such optical depths are simply not 
supported by the visual appearance of the Call K line in Fig|5|when 
compared to the hydrogen lines in Fig|2| 

As there are only two calcium lines, the uncertainty in the 
derived Call column densities are more difficult to estimate. The 
formal errors on each individual Call column density are almost 
certainly underestimates of the true uncertainties. The normalisa- 
tion of the continuum in these blue orders is likely to be the pri- 
mary source of error due to the difficulty of locating the contin- 
uum. The Call H&K lines are in subsequent echelle orders. We 
ran a Monte Carlo simulation assuming that the continuum of each 
of these orders may be inaccurate by up to ±5%. The standard 
deviation of the resulting probability distribution provides a sys- 
tematic uncertainty on all calculated Call column densities giving 
log Ncaii = 17.73 ± 0.25 m" 2 . A further check was made us- 
ing a template reference star taken during the observing run, GL 
472. This is a KOV dwarf observed as a slow rotating reference 
template for use in Doppler imaging. The narrow lines of this star 
make the continuum easier to fit. We therefore performed the same 
analysis, using instead the Call H&K continuum fits of GL 472. 
The resulting average is log Ncaii = 17.63 m~ , very similar 
to the result quoted above and within the systematic uncertainty 
determined above. 

We investigate the result of increasing the turbulent velocity 
from 0-20 km s~ x in Fig.|4] A vertical dashed line on the plot 
shows the maximum allowed turbulent velocity of 1 1 km s _1 . The 
resulting Call column density increases by a factor of five over the 
range 0-11 km/s from log Ncaii = 17.46 m~ 2 to log Ncaii = 
18.2 m~ 2 . This is discussed further in 35.21 

3.3 Weighing prominences 

From the previous sections we now have the column densities and 
optical depths for the four prominences in the n=2 level of hydro- 
gen and in the ground state of singly ionised calcium. In order to 
calculate the mass of each prominence we need to know the column 
density in the ground state of hydrogen. 

At the expected prominence temperatures of approximately 
10,000 K we assume that the calcium atoms are singly 
ionised. Then we can use solar abundance ratios to find the 
amount of neutral hydrogen. The ratio of calcium to hydro- 
gen in solar prominenc es is approximately log H g= 1 = 5.65 
( Gouttebroze & Heinzel 2002). We adopt this ratio to obtain col- 
umn densities in the ground state of hydrogen, displayed in Table 
□ 

We now require an estimate of the area of the promin ences. 
Here we use the same approximation as Collier Camero n et alJ 
( 1990). From Jy] 

we know that the the prominences are optically 
thick at K a (r ~ 22) line centre. Therefore when they transit the 
centre of the stellar disc, the fractional depth of the absorption pro- 
file gives us the fraction of the star obscured by the prominence. 



This is analogous to how the fractional depth of a planetary transit 
gives the ratio of the areas of the planet and star. We use this to cal- 
culate the absolute areas of the prominences, using the Speedy Mic 
radius of 1.06 ± 0.04 R© found in Paper I. These are shown in 
Table □ The average area is found to be 20% of the stellar disc. 
The masses of the prominences can now be calculated as simply 
the product of the hydrogen column density, the prominence area 
and the mass of a hydrogen atom: 

M = m H N 1 A 

The results are again shown in TableQ Prominences A, C and 
D are found to have very similar masses while prominence B has 
considerably higher mass. For a turbulent velocity of 5 km s _1 the 
average of the four prominences is 1.3xl0 14 kg. We note that the 
individual uncertainty on each prominence arising from systematic 
error (as discussed in 33. 2t exceeds the scatter in the prominence 
masses resulting in the range 0.7 — 2.3 xlO 14 kg. 

The Call column densities found in the last section, from con- 
sidering the full range of permitted turbulent velocities, provide an 
estimate of the minimum and maximum prominence masses. These 
are found to be 0.5 — 3.4 xlO 14 kg. We discuss these limits further 
in jB"2l 



4 PROMINENCES SEEN IN EMISSION 

Solar prominences are visible as bright emission features when they 
are seen off the limb of the Sun. In theory it should be possible to 
detect stellar prominences as bright emission features when they 
are off the limb of the star. The source function of a prominence 
will be stellar H a photons that have been scattered (see into our 
line-of-sight by the prominence, plus any additional source func- 
tion from the thermal properties of the cloud itself. The relative 
intensity we observe will be proportional to the size and optical 
thickness of the prominence in H a and inversely proportional to 
the square of the prominence height above the stellar surface. As 
Speedy Mic prominences are found at least two stellar radii above 
the stellar surface these features are likely to be very faint. It is also 
worth noting that, depending upon the geometry of the system, a 
prominence seen off the stellar disc need not transit the star as seen 
from our line of sight. 

4.1 Detecting loops of emission 

The exceptional signal-to-noise of the VLT H a timeseries allows 
us to examine the regions beyond the stellar usini limits for evi- 
dence of off-disc emission. To do this we created a mean profile, 
this time using the whole of the timeseries as there are no domi- 
nant strong features present outside of the stellar H Q profile. This 
was then subtracted from the entire timeseries to search for time- 
varying features. 

We first need to ascertain a conservative estimate for our sen- 
sitivity to weak off-disc features. The formal errors on each pixel, 
derived from Poisson statistics during the reduction, are underesti- 
mates as they take no account of systematic variations throughout 
the night. Thus a featureless section of the timeseries adjacent to, 
but unaffected by, the H Q line and also avoiding any strong iron 
lines was chosen. The standard deviation from the mean in this sec- 
tion which included all phases was a — 0.004, or 0.4% of the local 
continuum level. This gives us a guide to the significance of any of 
the features we identify. 

In the right-hand panels of Fig|6|we show the resulting H Q 
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Figure 5. The August 02 Call K (left) & H (right) timeseries are shown with contrast set so that black is at 0.85 and white is at 1.05 times the local continuum. 
An average profile has been subtracted (see text). Dashed vertical lines show the v sin i limits. Again, note that the Call H line is blended with the H e line. 
Bright surface plage regions are visible as well as weak off-disc emission. A stellar flare is present at phase fs 0.3. 



timeseries for both August 02 and 07, with contrast limits set to 
reveal small fluctuations above or below the continuum level. Fig[(>| 
shows a number of bright emission loops. In particular there is a 
very distinct loop of emission on August 02 to the red (positive ve- 
locity) side of the H a line, centred on <f) ~ 0.8. The emission peaks 
at 2.5% of the continuum level but weakens beyond the maximum 
elongation. The loop starts just after prominence A at <f> « 0.56 fin- 
ishes transiting the stellar disc. One immediate observation is that 
there is no corresponding emission loop on the blue side of the H a 
line. 

There are two other loops of emission both centred on <j> w 
0.15. The smaller loop shows very strong emission, peaking at 
3% of the continuum, but it cannot be separated from the H a line. 
This is probably a prominence very close to the stellar surface. The 
larger loop peaks at around 2% of the continuum level. Again it is 
apparent that there is no corresponding loop on the red side of the 
H Q line. The observations on August 07 are generally of a lower 
quality yet similar emission loops can be seen in Fig. [6] 



4.2 Understanding the emission 

In order to attempt to understand the nature of the relationship be- 
tween the on-disc absorption and off-disc emission features we can 
use what we have learnt from tracking the prominences in absorp- 
tion. For a detailed explanation of this matched filter technique we 
refer the reader to Paper I. We use the calculated phases and heights 
displayed in Table|2|to complete velocity sinusoids which we trace 
out in Fig|6| We do this for the three strongest prominences on the 
stellar disc for both nights. These are prominences A, B and C on 



Table 2. The results of the prominence tracking analysis are displayed with 
phases of meridian crossing and calculated heights. 



Prominence Phase Height ( ^S- ) 



A 


0.560 


2.99±0.07 


B 


0.378 


3.68±0.26 


C 


0.184 


3.72±0.20 


D 


0.406 


3.92±0.21 


E 


0.867 


2.35±0.07 


F 


0.258 


4.67±0.16 



August 02 and prominence D on August 07. We also include the 
next two strongest prominences on August 07, labelled E and F. 

In general there is good agreement between the extended si- 
nusoids and the observed loops of emission. For example promi- 
nences A and B on August 02 clearly account for the two strong 
loops of emission. Prominence C, however, does not seem to corre- 
late with any emission feature. A similar off-disc emission pattern 
is seen on August 07. Prominence D appears to account for most 
of the emission on both the blue and the red side of the timeseries. 
While the low lying prominence E appears to correlate with a loop 
of emission seen to the blue side of the line at <j> « 0.65. The high 
prominence F (like C) does not appear to be associated with any 
emission loop. 

Furthermore, the paths of prominences as estimated from 
tracking the absorption features appear to have systematically 
greater velocity amplitudes than the majority of the corresponding 
off-disc emission features. This can be explained in the context of 
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Figure 6. Timeseries to show loops of off-disc emission. The left panels show the fits to the transiting absorption features and have a grey-scale that runs from 
black at 0.85 to white at 1.05 times the local continuum level. Similarly the top right panel is set 0.98 to 1.02, while the bottom right panel has 0.99 to 1.01, 
in order to highlight the weak emission. The fits from the left panels are extended as velocity sinusoids (red) on the right panels. Dashed vertical lines (green) 
show the v sin i limits. 



the prominence optical depths as obtained in ij3| The average op- 
tical depth of these prominences was found to be r = 22.4. Thus 
as stellar photons travel through the prominence interior they are 
scattered into random directions, some into our line of sight. How- 
ever because of this scattering the beam loses intensity and so the 
number of photons reaching the highest parts of the prominence is 
reduced. 

When finding prominence heights from tracking the absorp- 



tion features we are naturally following the 'lack of photons'. Yet 
when we examine the off-disc emission we observe where these 
photons have gone. As such, the observed emission provides us 
with information on the radial extent (or thickness) of the promi- 
nences. Clearly from a visual examination of Fig|6|the prominences 
we observe do not extend to the stellar surface. We scale a Gaussian 
model to the most prominent loop of emission on August 02, fitting 
for the width of the feature. An example fit for a spectrum taken at 
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Figure 7. An example of a single H Q spectrum at maximum elongation of 
the strongest observed loop of emission seen at <f> »* 0.8 on August 02. The 
dashed lines show the stellar vs'mi limits. The solid line shows the Gaussian 
model we have fitted to the emission, which is found to peak at 2% above 
the continuum level. The FWHM is 260 km s _1 giving an estimate of the 
prominence thickness of 2.0 R*. 

maximum elongation of this emission loop is shown in Fig0 Both 
the amplitude and width of these emission features vary consider- 
ably along the emission loop. We find a range of full width at half 
maximum (FWHM) that corresponds to a prominence thickness of 
between 1.7 and 2.5 R». 



5 DISCUSSION 

5.1 Prominence stability 

In Paper I we investigated prominence stability on the timescale 
of days we found that many prominences could be re-observed 
on consecutive days, while a few were still recoverable some five 
nights later. In i|2.1l we extended the timebase by comparing Paper I 
observations with the current data taken more than two weeks later. 

The main result of this comparison is that the overall circum- 
stellar distribution of prominence material remains the same. While 
there is clearly evolution of the detailed structure of the individual 
prominences, the locations at which prominences are likely to be 
found are largely unchanged. In particular prominence A was the 
structure we could most easily identify between the two epochs 
of observations. We found that its phase of observation, height 
and strength were all remarkably similar. This does not necessar- 
ily mean that this feature has remained unchanged for two weeks 
(without observations of the intervening period it would be impos- 
sible to tell). What it does suggest however is that the magnetic 
structure within which this prominence was trapped is stable on 
these timescales. Within this structure, individual prominences may 
have formed and either been ejected or drained back to the surface 
within the two week gap when we do not have observations, but if 
this is the case, then another prominence must have reformed in the 
same place and under very much the same conditions. 

It is interesting that this prominence, which appears to be the 
most stable, i s also one of the l argest. This is also true of solar 
prominences fcchriiver & Zwaarj fcoOO ) . Small solar prominences 
which form close to active regions tend to be most unstable while 
the longest-lived are also typically the largest. In the solar case, 
the largest prominences are sheared out over their lifetimes by the 



action of differential rotation and they are carried poleward over a 
lifetime of months to form the polar crown. Their poleward mi- 
gration is halted when they encounter the open field of the po- 
lar cap. Whether something similar is happening on Speedy Mic 
is not clear. Certainly, the apparent lack of high-latitude promi- 
nences (discussed in £|5. 3> suggests that the prominences are con- 
strained to a range of latitudes, rather like their solar counterparts. 
This is in stark contrast to the distribution of spots on Speedy Mic 
which can be seen at all latitudes, including the pole lBarnesl2005l 
IWolter et all20ol . 

5.2 Prominence masses 

In J3|we used the ratio of the strengths of multiple Balmer lines 
to determine the column densities of prominences in the first ex- 
cited state of hydrogen t o be (log N2 = 18.61 ± 0.2 m ~ 2 ). This is 
in good agreement with lCollier Cameron et alJ ll99(t) who found 
(log N2 = 18.4 ± 0.2 m -2 ). Consequently similar optical depths 
are also derived. From direct measurements of the EWs of promi- 
nences in the Call H&K lines, and then use of solar prominence 
abundance ratios, we were able to estimate the column density in 
the ground state of hydrogen. This was found to be (log iVi = 
23.38 ± 0.25 m~ 2 ). Then using the projected areas of the promi- 
nences we estimated their masses to be 0.7 — 2.3 xlO 14 kg. 

The results outlined above were for a turbulent velocity of 5 
km s _1 . The turbulent velocity could be in the range 0-11 km s _1 
producing a larger prominence mass range of 0.5 — 3.4 xlO 14 kg. 
However there are reasons why we believe that the turbulent ve- 
locity probably lies towards the lower end of this range. Firstly we 
can examine the assumption that the smallest prominences are un- 
resolved. The smallest features for which we could determine the 
widths have fractional depths in H a that correspond to areas that 
cover at least 2% of the stellar disc. Assuming a circular cloud 
cross-section this gives a radius of 14 % that of the Speedy Mic 
radius. Hence in velocity space a width of 18 km/s would be ob- 
served. This in itself accounts for the observed widths. 

Another argument for small turbulent velocities comes from 
comparing the optical depths of prominences in the H Q and Call 
K lines. On Fig. |4]we include a second axis that shows the corre- 
sponding optical depth in the Call K line. For the range of turbulent 
velocities 0-11 km s" 1 optical depths of 8 < r < 40 are found. 
We also plot on Fig.|4]a horizontal dot-dash line to show the optical 
depth at the centre of the H a line at r = 25.7. If we now compare 
the Call K line in Fig.|5|with the that of the hydrogen lines in Fig. 
|2| the prominence absorption features in the Call K line appear to 
be less optically thick than the H a line. This includes the lack of 
an emission loop associated with prominence A. The Call K time- 
series is of sufficient S/N that we would expect to detect scattered 
photons from this loop if it had a similar optical depth to that of 
H a . We suggest then that the optical depth of the Call K line is 
less than the H„ line and thus the turbulent velocity should be less 
than 6 km s~ . This is therefore similar to the level of turbulence 
seen in solar prominences. Hence our adopted turbulent velocity of 
5 km s _1 seems valid in the previous sections. Therefore we arrive 
at a narrower range for prominence masses of (0.5 — 2.3 xlO 14 kg). 

We can now compare some of the physical properties of the 
Speedy Mic prominence system found in J3|with those of the only 
other similar analysis on a stellar prominen ce system. Usin g differ- 
ent techniques to this work. lCollier Cameron et all 179 901 derived 
parameters for another rapidly rotating star, AB Dor. They found 
prominence masses of 4 — 5xl0 14 kg, slightly greater than we 
find for Speedy Mic. Given the good agreement between this work 
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and lCollier Cameron et alJll99d) on the column density of the n=2 
level of hydrogen it seems unlikely that the above discrepancy on 
prominence masses represents a physical difference between AB 
Dor and Speedy Mic prominences. Rather it is probably just a dif- 
ference between the two techniques employed due to the quality of 
data available at the time. With the advent of new highly-efficient 
spectrographs and increased detector size and efficiency (particu- 
larly in the blue part of the spectrum) it is now possible to obtain 
simultaneous high-resolution spectra in the majority of the optical 
wavelength range. In particular we are now able to obtain simulta- 
neous observations of the H a line and both Call H& K lines. This 
was not possible at the time of the ICollier Cameron et alj Il990t) 
observations, they therefore had to rely on a geometrical argument 
to obtain estimates of the column densities. 

We can also compare Speedy Mic prominences with their 
solar namesakes. For example a gia nt eruptive solar prominence 
was observed by GoDalswamv & Hanaoka 1 1998) to have an initial 
mass before eruption of 6xl0 13 kg. This is comparable to the lower 
end of our mass estimate. However to have three such prominences 
on the Sun at any one time would be unusual indeed. 

As we have the column densities in the first two levels of hy- 
drogen we can estimate the excitation temperature of the promi- 
nences by assuming that they are in local thermodynamic equilib- 
rium (LTE). We note that this can not strictly be the case as else the 
10,000 K prominences would appear as emission rather than ab- 
sorption features whilst transiting the stellar disc. However given 
that we have shown H a to be just optically thick, we consequently 
know that Lyman alpha must also be very optically thick. As it 
is this transition that is largely responsible for detailed balance 
between the first two levels of hydrogen the LTE approximation 
should be reasonably accurate. 

Comparing the column densities of hydrogen in the ground 
state with that in the first excited state, the average shows that there 
are 10 4 ' 77 = 59, 000 times more atoms in the ground state than in 
the first excited state (assuming a turbulent velocity of 5 km s _1 ). 
Therefore the ratio of the LTE level populations is given by the 
Boltzmann equation: 




where gi and g2 are the statistical weights of the two levels, AE 
is the energy difference between them, k is the Boltzmann constant 
and T is the temperature sought. This results in an excitation tem- 
perature of 9569 K. This is very similar to the kinetic temperature 
of 10,000 K we used in section [J3] While we realise that this is a 
circular argument, the results of this analysis are not found to de- 
pend strongly on the thermal Doppler temperature used in creating 
the theoretical line profiles. For example if we start with a lower 
thermal temperature of 8000 K we obtain an excitation tempera- 
ture of 9463 K. 



5.3 Distribution of prominence material 

In !j4|we showed that some of the strongest absorption features (the 
largest prominences) could be observed as loops of emission. How- 
ever we found no emission counterparts for prominences C and 
F. This probably reflects the relationship between the size of the 
prominence and the height above the stellar surface. Here C and F 
cover slightly larger areas of the stellar disc than say prominence E. 
Yet because prominence E is very much closer to the stellar surface 
(7? = 2.35 R») than either C or F (R = 3.72 R» and 7? = 4.67 R, 



respectively), and the intensity drops as the square of the distance, 
it is seen in emission whilst prominences C and F are not. 

If prominences were evenly distributed in latitude around 
Speedy Mic then one would expect to observe loops of emission 
which can not be associated with prominences crossing the stellar 
disc. However this seems not to be the case, from Fig. [6] we can 
associate all the observed emission loops with transiting promi- 
nences. If large prominences, at similar heights to those observed, 
existed at latitudes of around 60° then they should be visible as 
emission loops. So this either suggests that Speedy Mic does not 
have such large high-latitude prominences, or that they must be 
extremely centrifugally flattened towards the equatorial plane. It 
should be noted though that if higher latitude prominences were 
smaller or indeed higher then they may exist without detectable 
emission. 

Similar conclusions have been reached before from the 
lack of prominence s observed on low-inclination rapid rotators. 
Ijeffries et all J 1994 found no prominences transiting the stellar 
disc around the rapidly rotating K dwarf BD+22°4409 which is ob- 
served at an inclination of 50 ± 10°. Unfortunately, as discussed in 
the the Speedy Mic inclination is relatively poorly constrained 
and restricts what we can determine. If the inclination is actually 
the 70° that both Doppler imaging teams found then even a promi- 
nence above 7? = 3 R*, located in the equatorial, plane would not 
be seen crossing the stellar disc at all. Indeed, these observations 
may favour a higher stellar inclination. 

Speedy Mic is not the first rapidly rotating young star to 
show rotationally modulated em ission outside of t he H a line. 
Two G dwarfs RX J1508.6-4423 iDonati et alJl200d) and AP149 
(Barnes et al. 2001) also show emission beyond their stellar v sin i 
limits which were attributed to a ring (or torus) of prominence 
material. These two stars are very different cases to Speedy Mic 
though. On neither star are prominences seen transiting the stellar 
disc, nor is there any evidence that some prominences eclipse the 
off-disc emission of others. This is probably as both stars are seen at 
low inclination angles (i = 35° and i = 30° respectively). Further- 
more the strength of the emission is very different. The amplitude 
of the emission in these systems is around 20% of the local con- 
tinuum as opposed to the 2-3% we observe on Speedy Mic. This is 
probably due in part to the increased ionising UV flux from these 
early-type G dwarfs. This coupled to the fact that the co-rotation 
radius on these stars is closer to the surface than in Speedy Mic 
means the prominences receive greater irradiation. This drives the 
H Q source function closer to the thermal equilibrium value than to 
the pure-scattering value. 

Our analysis of the Speedy Mic prominence system has con- 
sidered only individual prominences that we can track across the 
stellar disc. This is a simplification of the real distribution of promi- 
nence material. For example if we look back at the H a timeseries 
shown in Fig|2|we can see that prominences A and B are not to- 
tally separate entities. The H Q line is still in weak absorption at the 
phases between these two prominences. Yet if we compare this with 
the later Balmer lines (like H 7 in the panel below H a in Fig|2j the 
absorption weakens considerably. This indicates that this material 
has a lower column density than the material in what we refer to 
as prominences A and B. The large prominences that we have been 
studying therefore appear to be more like dense condensations in 
a 'semi-torus' of gas. This would also help to explain why some 
prominences do not show off-disc emission loops on both sides of 
the H Q line, as observed in J4] It is the intervening material be- 
tween some of the prominences that obscure us from observing this 
emission, i.e. the prominence system is 'self-eclipsing'. Therefore 
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there appear to be similarities between the Speedy Mic prominence 
system and that of the low inclination systems RXJ1508. 6-4423 
and AP149. The observed differences between these systems and 
Speedy Mic in many respects appear to be due to the different stel- 
lar inclinations. 



6 CONCLUSIONS 

The first analysis of the evolution of a stellar prominence system on 
a two to three week timescale has been performed. Whilst we see 
considerable evolution of most individual prominences the overall 
prominence pattern remains broadly similar. One half of the stel- 
lar rotation is still covered in stronger prominence signatures than 
the other. We have shown that at least one individual prominence 
supporting structure appears to be stable over two weeks. Any the- 
oretical models for prominence support will have to be able to re- 
produce such stability. 

We show that the prominence pattern is visible in the first eight 
lines of the hydrogen Balmer series and the Call H&K lines. The 
column density of hydrogen in the first excited state and the derived 
H Q optical depth ar e very similar to that of the prominences on 
AB Dor reported bv lCollier Cameron et alJ 1199(J> . This analysis, 
however, finds a lower total hydrogen column density from direct 
observations of Call. This leads us to find prominence masses of 
0.5 — 2.3 xlO 14 kg, intermediate between the masses of the largest 
solar prominences and those of AB Dor. We suspect that the differ- 
ence in the derived masses of the prominence systems on these two 
rapidly rotating stars is not physical but is probably a result of the 
different techniques used and the data quality available. 

We have identified loops of emission outside of the rotation- 
ally broadened H a profile. These are shown to be caused by the 
same prominences we see transiting the stellar disc in absorption. 
The emission loops confirm the heights found from tracking the 
absorption signatures as the visible emission is shown to appear 
at smaller radii in all cases. This observation can be explained by 
the optical thickness of the clouds in H Q . By combining these two 
methods of observing prominences we have been able to develop 
a more complete picture of the geometry and physical parameters 
of the prominence system. All the strong loops of emission we ob- 
serve can be associated with particular prominences seen transiting 
the stellar disc. This limits the number of large high-latitude promi- 
nences. Furthermore, combined with the high degree of self-eclipse 
observed in the Speedy Mic prominence system we suggest that it 
is indeed highly flattened, most probably near to the stellar equato- 
rial plane. This apparent flattening of the prominence system could 
result from the presence of a unidirectional field in the polar cap. If 
there is a strong dipolar component in the stellar field then promi- 
nences will form in the equatorial plane l Mclvor et al. 2003). From 
Doppler imaging of the Speedy Mic surface 16*^16X12003) found 
that there was no single uniform polar spot. So it may be the case 
that the stellar field has a high degree of complexity and promi- 
nences that form at high latitudes are centrifugally driven into the 
equatorial plane. 

It is clear that the Speedy Mic prominence system is complex, 
with the distribution of co-rotating prominence material not re- 
stricted to the dense clumps that we have identified as prominences. 
We note that there seems to be interesting parallels between these 
observations and those of the low inclination G dwarfs RXJ1508.6- 
4423 and AP149. In the future it may be informative to develop 
a model where we could change the physical parameters and the 



viewing angle of such prominence systems. The results could then 
be compared directly with these observations. 
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